Singular length dependence of critical current in SNS bridges 



o 
o 

Or 
(N 



Alex Levchenko, 1 Alex Kamenev, 1 and Leonid Glazman 1 ' 2 

1 Department of Physics, University of Minnesota, Minneapolis, MN 55455, USA 
^W.I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

(Dated: October 25, 2006) 

We examine dependence of the critical Josephson current on the length L of the normal bridge N 
between two bulk superconductors. This dependence turns out to be nonanalytic at small L. The 
nonanalyticity originates from the contribution of extended quasiparticle states with energies well 
above the superconducting gap. This should be contrasted with the more familiar contribution to 
the Josephson current coming from Andreev bound states localized in the normal region at energies 
below the gap. 

PACS numbers: 74.50.+r, 74.45.+C 
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Owing to the rapid progress in nanofabrica- 
tior*ii£i^, there is a renewed interest in various 
aspects of the physics of superconductor/normal- 
metal/superconductor (SNS) structures^. The purpose 
of this paper is to point out a subtle contribution to 
the Josephson current, flowing across the normal region. 
Specifically, we address the current carried by the high 
energy, e > A, extended states of the system. It leads 
to a nonanalytic behavior of the critical current as a 
function of the length of the normal region, L, at small 
L, 

It is frequently stated^ that as the result of Andreev 
reflection^ the Josephson current is carried exclusively 
by the Andreev bound states with the energies e < A, 
localized in the normal region. Indeed, the energies of 
such bound states are sensitive to the phase difference, 
</>, between the superconductors. This leads to the 4>- 
dependent free energy and thus to the supercurrent. The 
validity of this point of view is well established in the 
two limiting cases of lon g 9 : 10 : 11 junction L ^> £, and very 
short£i£i£ one, L/£ — > 0, where £ is the coherence length 
of the superconductor. Thus it came as a surprise for us 
that in between these two extremes, the physics appears 
to be more complicated. We find that the Josephson 
current is shared between the localized Andreev states 
and the extended above-gap states. The dependence of 
the critical current on L at short lengths, L <C £, is not 
analytic and comes from the contribution of the extended 
states. 

The phase sensitivity of the energy levels originates 
from the trajectories, which are reflected at least once 
from both NS interfaces. Such trajectories, require prop- 
agation time longer than the diffusion time across the 
normal region L 2 / D, where D is the diffusion constant. 
As a result, only states with energy smaller than the 
Thouless energy, e < E^h = hD/L 2 , may exhibit sensi- 
tivity to the phase difference. In the long junction limit, 
the Thouless energy is small compared to the supercon- 
ducting gap, i?Th "C A (indeed, notice that A = HD/^ 2 ). 
Thus all the current-carrying states are at the energy 
e < Erh "C A. Since there are no states in the supercon- 
ducting leads at such energies, the Andreev bound states 
of the normal region arc indeed solely responsible for the 



Josephson current. 

The situation is more complicated for shorter junc- 
tions, such that i/Th > A. One may expect that the 
extended states in the energy range A < e < Et\i 
contribute to the supercurrent. Yet, in the only case 
treated analytically so far, namely in the limit L/£ —» 
(or, equivalcntly, i?xh - ► oo) the entire current is still 
given by the Andreev states^ with < e < A. In 
this paper, we show that this result is an artefact of the 
limit i?Th - ► oo rather than a general principle. Specif- 
ically, the spectral density of the Josephson current ap- 
pears to be J(e) ~ A 2 /(ETh£) in the energy interval 
A < £ < -Erh- Upon integration over energy, it results 
in the contribution to the critical current of the form 
SI C ~ (A 2 /£/rh) ln(A/£/rh), nonanalytic in the limit of 
the short normal region. 

Summarizing the existing knowledge, one may write 
for the critical current of an SNS structure 



I c = — 1(k), 
e 



A 



3 



Th 



(1) 



where g is the conductance of normal region and T(k) is a 
universal scaling function of the junction's dimcnsionlcss 
length k = (L/t;) 2 . In the limit of long junction, n ^> 1, 
this function is given byiL T(k) ps 10.82 On the 

other hand, for shorter junctions, k <C 1, it is 



I(k) = 1(0) - an\n(b/K) + 0(k 2 ) , 



(2) 



where J(0) = 2.08, see, e. g., Ref. 0, and a = 0.31, 
while b — 2.84. The first term on the right-hand side 
(r.h.s.) represents the current carried by the Andreev 
states and was discussed by Kulik and Omelyanchuk 8 
and Beenakker— . The second nonanalytic term originates 
from the above-gap extended states. This term is the 
main result of the present paper. In the rest of the paper, 
we derive it and compare our analytical result with the 
existing^ numerical data. 

To approach the problem analytically, we employ the 
Usadel equation^ for the matrix Green function G(e, x) 
of the quasi-lD disordered normal region, 



Dd x (Gd x G)+i£[a 3 ,G}=0 : 



(3) 
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where the Green function is subject to the constraint 
G 2 = 1. At the two NS interfaces, the Green functions 
G(e, ±L/2) are given by those of the BCS superconduc- 
tors with the order parameter A, maintained at the phase 
difference <p. The suppression of the order parameter 
within the bulk superconducting leads may be safely dis- 
regarded as long as the transverse directions of the nor- 
mal region are less than the coherence length, £. The 
same assumption justifies keeping only the lowest trans- 
verse mode and thus the use of the ID form of the Usadel 
equation. This should be contrasted to the case of the 
long junction L 3> £, where fully self-consistent calcula- 
tion is required and the approximation with the constant 
order parameter within the normal region is not appro- 
priate. 

The Green function may be parametrized by the two 
complex angles 6(e,x) and p(e,x) as 

G = sin 9 cos tp ai + sin 9 sin (p &2 + cos 9 a 3 . (4) 

Employing this paramctrization and introducing the 
rescaled coordinate x/L — > x, we rewrite the Usadel 
equation as 



dl9 + uj 2 sin 9 = J 2 cos 9 (sin 9) 3 , 



(5) 



where ui 2 = 2ie/£ , Th, and the spectral current density 
J(e) is given by 

J = sin 2 9d x p . (6) 

One needs to solve Eq. ([5]) with the boundary conditions 

tan0(e,±l/2) = iA/e. (7) 

The spectral current density may then be determined 
from the condition of having the fixed phase difference: 
(p(e, 1/2) — tp(e, —1/2) = </>. Finally, the Josephson cur- 
rent is found as 



00 



(8) 



For the sake of illustration, we shall first execute this pro- 
gram in the short junction limit, L ~ uj — ► 0. Equation 
([5]) with uj = may be easily integrated, resulting in 



cos 6(e, x) = cos 9q cos 



Jx 
sin 9r 



(9) 



where the integration constant is 60(e) = 9(e,0). Inte- 
grating then Eq. ^ and employing the boundary condi- 
tion for the phase tp, one obtains 



tan(<£/2) 



1 



sin 9q 



tan 



J 



2 sin #0 



(10) 



This last equation along with Eq. ((9]) taken at the NS 
interface, x = 1/2, constitute the system of the two alge- 
braic equations for the two unknown quantities: J(e) and 



60(e). Such an algebraic problem may be easily solved, 
resulting in the following expression for the imaginary 
part of the spectral current: 



Im J = 



ttA cos(</>/2) 
^Je 2 - A 2 cos 2 (0/2) 



(11) 



for Acos(</>/2) < e < A, and Im J = otherwise. The 
fact that Im J vanishes for e > A is an artefact of the ap- 
proximation Etii — > 00, already mentioned in the begin- 
ning of this paper. Equation (jlip is in perfect agreement 
with the result of Beenakker— based on the consideration 
of the Andreev bound states in the diffusive normal re- 
gion. One concludes that the Josephson current in the 
limit L — ► is entirely given by the Andreev states resid- 
ing inside the superconducting gap. Employing Eqs. ([8]) 
and (jTTJ) one arrives at the well-known result^ for the 
zero-temperature Josephson current of the short SNS 
junction, 

I(4>) = I cos(0/2) arctanh[ sin(0/2)] . (12) 

with I = irgA/e. At ma x = 1.97, this function reaches 
its maximum value I c = I(4>max) = 2.08<?A/e, which 
defines the coefficient 1(0) = 2.08 in Eq. ©. 

Having established the limit of an extremely short 
junction, we turn now to our main subject: the finite 
length correction to the critical current. As was men- 
tioned in the introductory section and will be proven 
below, the largest correction originates from the para- 
metrically wide range of energies well above the super- 
conducting gap: A < e < i?Th- One may notice that 
in this energy interval, the solution for 9(e) must be of 
the order of A/e 1 , allowing for small-0 expansion in 
Eq. ©, 



0. 



lu 2 9 = J 2 9~ 3 



(13) 



Within the same approximation, the boundary condi- 
tions © read 6(e,±l/2) = iA/e, while Eq. © for the 
spectral current density takes the form J = 9 2 d x (p. 
Equation Q13p may be solved exactly, leading to 

/v2/ \ n-2 ?/ x J 2 sm 2 (ojx) . . 

9 2 (e,x) = 9 2 cos 2 (ujx) + - t Y-i-, (14) 

"0 w 

where 60(e) = 9(e,0) is the integration constant, similar 
to the one in Eq. @. Substituting the solution Eq. (THl) 
into the equation for the spectral current and integrating 
over the coordinate, one finds 



tan((/>/2) 



J tan(w/2) 



(15) 



Notice that in the limit u> — ► Eqs. (Ti"4"|) and (|15p reduce 
to the small #0 and small J limit of Eqs. ^ and (|10[) . 
respectively, as they should. 

Taking Eq. (TIT} at the NS interface, x = 1/2, and 
employing the boundary conditions, one obtains the al- 
gebraic relation between yet unknown quantities 9o and 
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J. The second relation between these two quantities is 
provided by Eq. (fT5| . Resolving these two algebraic re- 
lations, one finds for the spectral current density 

A 2 oj 

J = — — sin<p . (16) 

e z sin to 

In the short junction limit, to — > 0, this expression is 
real, which is in agreement with our previous finding 
that Im J = for e > A. How ever, for a finite-length 
junction (recalling that to = ^J2\e/Erh) one finds the 
non-zero imaginary part of the spectral current: Im J w 
— A 2 /(3£ , Th^) in the energy range A < e < £rh' For 
larger energies, e > -Erin Eq. (fTo| predicts exponentially 
decaying Im J oc cxp(— yJe/E^). Indeed, as was dis- 
cussed in the introductory part, the energies above the 
Thouless energy should not contribute to the supercur- 
rent. Note that the tail of the spectral current found in 
the Eq. (JT5J) at the energies above the gap is in perfect 
agreement with the results of the numercial studies^. In- 
tegrating Eq. (flB]) over the energy according to Eq. j8]), 
one finds the finite-size correction to the T = Joseph- 
son current 



14 



SI(</>) = -I 



= -tJ- 



V3tt 



ln(fe/K) sin^> . 



(17) 



This correction modifies the value of the critical current, 
5I C = <5/(0 max ), with max quoted after Eq. (|12|). This 
way we obtain the second term in the r.h.s. of Eq.© 
with the coefficient a = (sin0 max )/3 = 0.31. 

We comment here on the additional approximation 
made while deriving Eq. (|17p . Our solution strategy was 
based on the perturbation theory for the Usadel equa- 
tion over L/£ -C 1 exploring smallness on the phase 
6(x,e) in the energy range A < e < -Erh- But we have 
used rigid boundary conditions ([7|) ignoring corrections 
to them coming from the finitcness of the normal region 
link. Although these corrections exist, they do not bring 
any new contributions non- analytical in L/£ to Eq. (| 1 7|) . 

The crossover between the limits of short and long 
junctions was recently studied numerically by Dubos et 
ai.—. In Fig. [TJ we use the data of Ref. LU] to plot the 
scaling function 2(n). At small «, we expect the linear 
dependence of [X ~ T(0)]/k on ln/c, see Eq. ©. The 
numerical data^i indeed agree with this expectation, as 
shown in the inset to Fig. [TJ 

All the above consideration was based on the assump- 
tion that the interfaces between normal metal and su- 
perconductor are perfectly transparent. In most experi- 
mental configurations, SNS junctions always contain the 
potential barriers, an insulating layer I, at the SN inter- 
faces. Thus it is of interest to investigate the influence of 
the extended quasiparticlc energy levels on the Joseph- 
son current in the realistic SINIS configuration. In what 
follows, we will constraint ourselves to the case of sym- 
metric SINIS junction with tunnel barrier conductance 
<?t being much smaller then g. 

The physics of the Josephson effect in the SINIS junc- 
tion is controlled by the dimensionless interface param- 
eter 7off = (L I '£) 2 (g I gr) ■ In the so-called incoherent 




FIG. 1: (Color online) Numerical results of Ref.— for the 
scaling function Inset: to facilitate the comparison of 

numerical results (Ref.fTlT) (data points connected by a dashed 
line) with analytical theory [solid straight line corresponding 
to Eq. the data are replotted in co-ordinates [I— 1(0)] /an 
vs \n(n/b) for the region k < 1. 



regime 7 c fr 3> 1, the SINIS junction may be viewed as 
the two tunnel junctions connected in series with the 
current-phase relation in the form I(4>) = I c sin (0/2)^. 
The most interesting case is the coherent regime j c s <C 1 
which we consider below. 

To calculate the supercurrent for the SINIS junction, 
we start from the point contact limit and follow the same 
steps (p3T)- (fTU)) as for SNS case but with one significant 
difference, namely the boundary conditions for Usadel 
equation. Instead of rigid boundary condition ([7]), we 
apply conditions appropriate for this case^ 



2 9T ■ a ■ a ■ 

sin v g sin t> s sin 

9 2 



(18a) 



2ffr 

Ox0\ x= l = 

9 



sin 9s cos Ob cos ■ 



Vb 



sin 6b cos 6 c 



(18b) 

where 6b(s) — 6(e,±l/2), ips = '■P^A/'^) and cos#s = 
e/\Je 2 — A 2 . As long as the tunnel conductance is small 
9t g the phase function tp(e,x) changes abruptly at 
the tunnel barriers from its values ±</>/2 inside the su- 
perconductors to iipB at the boundaries of the junction 
and remains small inside the normal region \ip(e, t)| < 1. 
These observations implies that in the boundary condi- 
tions (|18|) one can safely put tps = 0. As follows from the 
Eq. (fTB|) the spectral current J is small as ~ <?t/<? <C 1. 
Therefore, according to the solution ^ of the Usadel 
equation the #(e, x) is essentially coordinate independent 
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d x 9(x,e) — 0, provided i?Th ^> Ag/g-r- In this case 
the boundary conditions (fT5|) represent the closed system 
of algebraic equations for unknown functions 9b and </, 
which can be solved in terms of <f> and 9s- As a result the 
imaginary part of the spectral current takes the form 

W = gl A2sin ^ (19) 

g VA 2 - e 2 v / e 2 - A 2 cos 2 (0/2) 

Combining this with Eq.® at zero temperature, we find 
the Josephson current-phase relation in the for m 15 16 

I{4>) = I {2gr/irg)K[axi(<f>/2)] sin0, (20) 

where K(x) is the complete elliptic integral of the first 
kind. Let us compare this result with Eg . (fT2"|) . First of 
all, we observe that the presence of the tunnel barriers 
changes the current-phase relation, but preserves its es- 
sential properties, for example nonanalyticity at phase 
<t> = 7r. Secondly, the amplitude of the critical current is 
suppressed by the small parameter gr/g <C 1 compared 
to its value, Eq. (fT2|) . Similar changes occur with the 
Josephson current correction SI: we get modified phase- 
4> dependence, the amplitude of the current correction is 
suppressed due to the tunnel barriers, but the correction 
itself remains nonanalytical in The estimate shows 
that the amplitude of the Josephson current correction 
in the symmetric SINIS junction is 

<5J c oc -I {g T /g)K\n{l/K). (21) 

Let us discuss now some characteristic features of the 
above-gap contribution to the Josephson current found 
in this paper: (i) The power-law tail of the spectral 
current density, Im J ~ 1/e, results in the nonanalytic 
length dependence of the critical current, 8I C ~ L 2 \nL. 
(ii) The supercurrent carried by the high-energy ex- 
tended states is negative, i.e., it flows in the direction 
opposite to the current produced by the Andreev bound 
states, (iii) Unlike the dependence on length, the phase 
dependence of the found contribution to Josephson cur- 
rent is not singular, SI(cf>) ~ sin0. This should be con- 
trasted with the contribution coming from the Andreev 
bound states, Eq. (|12[) . The latter contains (at any L) a 
nonanalytic phase dependence, /(</>) ~ (tt — 4>) In |tt — <f>\, 



at <f) ~ 7T. The origin of such a difference between the two 
contributions is in the evolution of the Andreev states 
with the phase <j> : unlike above-gap states, the lowest of 
the bound states depends on the phase in a peculiar way, 
"touching" the Fermi level, e = 0, at <fi = tt. (iv) The 
occupation of the states contributing to the L 2 In L de- 
pendence of the Josephson current hardly changes with 
temperature, as long as T < A(T). As a result, the 
temperature dependence of the second term in the r.h.s. 
of Eq. @ comes solely from the temperature-dependent 
superconducting gap A(T). (v) At a temperature close 
to the critical one: T c - T < T c < E Th , the short SNS 
bridge is in the regime where A(T) <C T. Performing 
energy integration according to Eq. one finds for the 
Josephson current: 

(22) 

where the first term in brackets originates from the An- 
dreev states, while the second one originates from the 
high energy extended states. The nonanalytic length de- 
pendence exists in this case as well. 

Although our calculations were carried out for the case 
of the diffusive normal region, it is clear that the effect 
is rather generic. In particular, with the proper redefini- 
tion of the Thouless energy (and possibly with different 
numerical constants), Eqs. (Q]) and @ should hold for 
ballistic SNS structures as welli&. Moreover, nonanalyt- 
ical correction survives even in the case of SINIS junc- 
tion but with suppressed amplitude due to small tunnel 
barriers conductance. Also the critical current of super- 
conducting weak links of the type SS'S should contain 
the nonanalytic term, originating from the high-energy 
states, and numerical studies^ only support this conclu- 
sion. 
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